Deconfinement phase transition in a two-dimensional model of interacting 2x2 

plaquettes 
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A two-dimensional model of interacting plaquettes is studied by means of the real space renormal- 
ization group approach. Interactions between the plaquettes are mediated solely by spin excitations 
on the plaquettes. Depending on the plaquette-plaquette coupling J, we find two regimes: 

"confinement" Jc < J < 1, where the singlet ground state forms an infinite ("confined") cluster 
in the thermodynamical limit. Here the singlet-triplet gap vanishes, which is the signature for long 
range spin-spin correlators. 

"deconfinement" < J < Jc, where the singlet ground state "deconfines" - i.e. factorizes - into 
finite n-clusters of size 2" x 2", with n < nc{J). Here the singlet-triplet gap is finite. 

The critical value turns out to be Jc = 0.4822... 

PACS numbers; 71.10.Fd,71.27.+a,75.10.-b, 75.10.Jm 



I. INTRODUCTION 

We will discuss in this paper the 2D Hamiltonian 

H = Ho + J-Hj (1.1) 

where Hq is given by isolated plaquettes occupied with 
spin-1/2 states and Hj describes nearest neighbour inter- 
actions of these plaquettes as shown in Fig. [1] For J = 1, 
the model (|l.ip reduces to the well known 2d antiferro- 
magnetic Heisenberg model. The modified Hamiltonian 
(II. ip has been proposed in studies of structural instabil- 
ities of two-dimensional systems. 

The singlet-triplet gap ("spin gap") has been studied 
by various methods: 

• Nonlinear a model as a low energy effective theory 
(0), i) 

• modified spin wave theory (Q), (jj) 

• cluster expansion up to fourth order starting from 
J — 0, i.e. isolated plaquettes (0) 

In the spin-1/2 case, the model (|l.ll) is expected (0) to 
have a quantum phase transition at a critical value Jc^, 
which is signalled by a vanishing singlet-triplet gap for 

J > Jc- 

The results of these works are: 

• no phase transition for the spin 1/2 case in the 
nonlinear a model (8) for any J < 1. 

• Jc = 0.118 for the modified spin wave theory [ver- 
sus Jc = 0.112 in linear spin wave theory (0)]. 

• Jc = 0.54 for the cluster expansion. 



• A value Jc = 0.555 was obtained by means of Ising 
series expansions (0)- 

• Jc has been determined by means of the CORE 
method (contractor renormalization expansion) 
first by Capponi et al. (0) and recently by Albu- 
querque et al. (fioh. These results are somewhat 
lower Jc = 0.548. 



• Recent Monte Carlo simulations ( TllfT^ yield val- 
ues close to Jc = 0.549. 

The nonlinear a model approach yields different results 
for different cut-off schemes. In [Kawakami et al. (8)] 
no phase transition was found for S = 1/2, whereas in 
[Takano et al. ([H)] a critical Jc (0.2 < Jc < 0.25) was 
obtained. 

The phase transition of the magnetic system has been 
discussed in its correspondence to a superfluid-insulator 
transition of the boson model ([13) . 

The authors of ref. (fToh start from singlet {S = 0) and 
triplet {S — 1) plaquette states. Excited states \S,m), 
S = 0,1,2, m = —S,..,S are absolutely necessary to 
generate interactions, since singlets alone cannot interact 
due to total spin conservation. 

In a recent paper (jT§), we have studied how interac- 
tions on the 4 plaquette compound - depicted in Fig. [1] 
- are created by single plaquette excitations. The con- 



^ The critical value Jc is related to the parameter 7c introduced 
inQby Jc = (l-7c)/(l+7c) 



FIG. 1: 4-plaquette system with 2" x 2" = 4 x 4 sites - here 
n = 2; the single plaquette on the right shows the enumeration 
of plaquette sites. 
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servation of total spin at each interaction point is imple- 
mented by means of the Wigner-Eckart Theorem for the 
transition matrix elements 



{S'i,m[\Sq{x)\Si,mi) 
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M{Six,Si). 



(1.2) 



They can be expressed in terms of a Clebsch-Gordan co- 
efficient and one reduced matrix element M{S[^x^ Si). 
The latter only depends on the initial and final plaquette 
spin Si, S[ and the triplet operator Sq{x) at site x. The 
phase Vq {v+ = — 1, ?;o=f- = l) results from the trans- 
formation properties of the spin operator Sq{x) under 
the group SU(2). The interaction between neighbouring 
plaquettes can be traced back to the product of reduced 
matrix elements at sites x and y (Fig. ^ 

X Y 



EHZH - EDIH 



FIG. 2: Interaction of neighbouring plaquettes. 

M{S'i,Si;S',,Sr) = Y.M{S'i,x,Si)M{S',,y,Sr) 

(1.3) 

Quintuplet excitations [S = 2) have not been consid- 
ered in ref. (fiol'). We want to stress here that triplet- 
quintuplet transitions are large - comparable with singlet- 
triplet transitions. It is shown in (flsh that the gaps 
(singlet-triplet and triplet-quintuplet) decrease in the 
renormalization process. We will see in this paper that 
the inclusion of quintuplet excitations will move the crit- 
ical value Jc substantially to a lower value. 

The paper is organized as follows: 

In Section[n]we summarize the details of the real space 
renormalization group approach in 2D models. 

In Section [ml we evaluate the renormalization group 
flow for various couplings and gaps. 

In Section IIVI we discuss the deconfinement of the 
ground state wavefunction for J < Jc- 

In Section |V] we present a finite size analysis of the 
singlet-triplet gap in both regimes: Jc < J < ^ (con- 
fined), < J < Jc (deconfined). 

Section |Vl] is devoted to the staggered magnetization. 



II. 



REAL SPACE RENORMALIZATION GROUP 
IN 2D MODELS. 



(2) 



In fill) we first studied the interaction matrices Aq 
of the four plaquette system (Fig. [T]) in the sectors with 
total spin S. The elements of the interaction matrices 
are fixed on one hand by the Clebsch-Gordan coefficients, 
which arise in the construction of eigenstates with total 



spin S (on the 4-plaquette system) and the evaluation 
of the Wigner-Eckart Theorem (|1.2p for the transition 

(2'] 

matrix elements. On the other hand only depends 
on the following couplings 



7 = - ■M(21;10) 

a 

13 ^ - •M(ll;ll) 

a 

e = - •M(22;ll) 
a 



and gaps 



P = 



E^ — En 



E2 + £'0 ~ 2Ei 



(2.1) 
(2.2) 
(2.3) 

(2.4) 
(2.5) 



We have factored out from the couplings (|2.ip - (|2.3p the 
"fundamental" interaction 



a = M(1,0;1,0) 



(2.6) 



which is induced by the singlet-triplet transitions on the 
plaquette. 

In Appendix A of ref. (fisl ) one can find the explicit form 
of the interaction matrices A5 (for J = l!) 5* = 0, 1,2 
under the premise that on the four plaquettes only ro- 
tational symmetric configurations with singlets, triplets 
and at most one quintuplet contribute. In this case the 
dimensions ds of the interaction matrices As turn out to 
be 

(do,rfi,d2) = (7,9,14). (2.7) 

The factor J in p.ip is taken into account in the in- 
teraction matrices A5(4, 4,7,/?, e) by a rescaling of the 



normalized gaps (|2.4p . (j2.5p whereas the couplings (|2.ip - 
()2.3p remain unchanged. 

Having constructed in this way the interaction matri- 
ces A^"-* on an n = 2 cluster (2" x 2") from the ground 
states on an n = 1 cluster, we turned to the question, 
whether it is possible in general to construct the inter- 
action matrix Ag"^^\ S = 0,1,2 from the correspond- 
ing quantities of a n x n cluster. This is indeed possi- 
ble under the assumption that the low energy states on 
the (n -I- l)-cluster can be built up again solely from sin- 
glet, triplet, quintuplet ground states on n-clusters. The 
n-dependcnce only appears in a renormalization of the 
couphngs (|2.ip - (|2.3p and energy differences (|2.5p . (|2.5p 

Here, we refer to dH) [eqns. (6.1)-(6.4); (6.5), (6. 6)] 
for the used formulas of the renormalization of the cou- 
plings and recursion formulas for the scaled energy dif- 
ferences. Note, that J does not appear in the first group 
of equations, whereas the remaining two (scaled energy 
differences) are linear in J. 

Each step n — > n -I- 1 in the renormalization procedure 
demands the diagonalization of the interaction matrices 
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a1?+'\ 5 = 0,1,2: 



= (2.8) 

^("+l)|^(n+l)^ ^ ^(n+l)|^(„+l)^ (2.9) 
^(n+l)|^(„+l)^ = (2.10) 



The eigenstates |cr("+i'), |r("+i)), 1^^"+^)) with the 
largest eigenvalues (t*^"+^-', r("+^\ enter in the 

quantities 

(a, 6) , G("+i)(a,6) (a, 6) = (1, 0), (2, 1) 
Fi"+i)(a,a) , Ff+iHa,a) (a, a) = (1, 1), (2, 2) 

according to the bilinear forms: 

= (2.11) 

G("+i)(a,6) = Ec["+'^G,,fc(a,5)rf+^) (2.12) 

i^("+i)(a,a) = ^(r("+i))V,,(a,a) (2.13) 

A; 

i^("+^Ha,a) = Y.{t^'^yFiAa,a). (2.14) 



The contraction /j^i(l,0), etc. are independent of n and 
listed in Appendix B of paper (fill). 



III. NUMERICAL EVALUATION OF THE 
RENORMALIZATION GROUP FLOW. 

We now turn to the numerical evaluation of the recur- 
sion formula of the couplings [eqns. (6.1) — (6.4) in ([Tsh] 
and gaps [eqns. (6.5), (6.6) in (15)] in order to study 
the n-dependence (i.e. finite size 2"+^ x 2"+-'^) and J- 
dependence. We start with the singlet-triplet gap p("+^\ 
which yields the signature for long range order: From Fig. 
[3] we see, that there are two different regimes: 

a) Jc < J < 1 

Here the singlet-triplet gap approaches zero with 
increasing system size. For J = 1 we are close 
to zero already on small systems for no = 3. For 
decreasing J, nQ{J) increases and seems to diverge 
for J ^ Jc- 

b) Below this critical value (J < Jc) the singlet-triplet 
gap does not converge to zero anymore. Note 
also that there is a change in the curvature of 
with n, which is for large n convex if Jc < J but 
concave if Jc > J. This allows for a very precise 
determination of 



Jc = 0.4822.. 



(3.1) 
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FIG. 3: The scaled singlet-triplet gaps p("+^' [(6.5) in Q^)] 
as function of n and J. 



Let us next turn to the coupling ratio r-r- [(6.1) in 

J— K • . . 2a(") ^\ ' 

(|l5f)]. As function of n this quantity has a maximum, 

which travels to larger values of n, if J is lowered (Fig. 
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FIG. 4: The coupling ratio [(6-1) in HI)] as function 

of n and J. 

For Jc < J < 1 all curves approach a common limit 
for large n: 



a 



(n+l) 



2a(") 



0.52. 



(3.2) 



For J < Jc we observe a monotonic decrease to a limiting 
value, different from 



2a(") 



= 0.25. 



(3.3) 



The n-dependence of the coupling ^^("+1) [(6.2) in (15)] 
is shown in Fig. [5l 

For Jc < J < 1 all curves approach a common limit 
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(«+i) 



1.0849. 



(3.4) 
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FIG. 5: The coupling 7("+i) [(6.2) in yji)] as function of n 
and J. 



for large n, whereas we observe a monotonic increase with 
n for J < Jr and a common limit 
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(n+l) 



1.3411. 



(3.5) 



We only want to mention that the "diagonal" couplings 
[(6.3) and (6.4) in (15)] , which do not change the pla- 
quette spins, die out after a few steps. 



IV. DECONFINEMENT OF THE GROUND 
STATE WAVEFUNCTION. 

It was pointed out in Section |TT1 that the renormaliza- 
tion group procedure demands in each step n ^ n + 1 the 
diagonalization (|2.8p - (|2.10p of the interaction matrices. 
We only keep those eigenvectors (|o'^"+^'')) with largest 
eigenvalue (cr^""*"^^). 

We want to discuss now the physical meaning of the 
eigenvector components: 

af+i) = (z,0;n+l|a("+i)) i = l,..,7 (4.1) 

in the orthonormal basis |i, 0; n + 1) in the singlet sector 
- defined in Table II of ref. dH). E.g. has 

to be interpreted as the probability to find in the singlet 
ground state the four plaquette configuration 








|1,0) = 

with four noninteracting singlets. 



If 



(n+l) 



(4.2) 



(4.3) 



the singlet ground state |cr'^"+^)) "deconfines" to the con- 
figuration 1 1,0) with four noninteracting singlets. 

In Fig. [51 we show (a^^^) as function of n - 
system size (2"+i x 2"+^). 



the 



FIG. 6: The probability {^^-i^^^^^ to find the configuration 
(|4.2|l with four noninteracting singlets in the ground state. 



For Jc < J < 1, 



(n+l) 



(n+l) 



converges to 



0.110236. 



(4.4) 



For J < Jc the deconfinement limit (|4.3p is practically 
reached at a finite value n = nc{J), cf. Fig. [3 

Tic (J) decreases with J < Jc and defines the largest 
cluster size (2"'=('^) x 2"=^'')) which is still confined. 

In summary we can say: In the two-dimensional system 
of interacting plaquettes defined in (jl.ip we observed two 
phases. 

• For Jc < J < 1, there is a confined phase, where 
the ground state does not factorize into finite n- 
clusters (2" x 2") but forms one infinite cluster in 
the thermodynamical limit n ^ oo. The vanish- 
ing of the singlet-triplet gap is the characteristic 
signature of this phase. 

• In the deconfined phase < J < Jc the ground 
state factorizes into finite clusters n < nc{J) where 
nc{J) defines the maximal size of clusters and is 
shown in Fig. [71 



V. FINITE-SIZE ANALYSIS OF THE 
SINGLET- TRIPLET GAP. 

In the confined regime Jc < J < 1 the singlet-triplet 
gap 



vanishes with an exponent 



1^1 



log(l -I- x) 
log 4 



(5.1) 



(5.2) 
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common limit 



FIG. 7: nc{J): as a criterion we show for each chosen J- value 
that number of iterations (n) where ^fj"^^'^ for the first 
time exceeds 0.999. 



which can be determined from the first derivative 
d(r — cr) 9(t — cr) 9(r — cr) dn 



dp 



dp 



dn 



dp 



(5.3) 



of the largest eigenvalues cr, r of the interaction matrices 
As{p, K, ...) in the singlet (5 — 0) and triplet (5 = 1) 
sector: 



i9(r - a) 



Ok 



9Ai 9Ao 

(r|^|r)-(a| — |.) 

-4(Tf - (7?) - 2(t| + t| + r| + r|) 
+2((7| + a2) (5.4) 



9k 



dn 

2 2 2 2 /c c;^ 

(EH) - which is the 



The n- and J-dependence of — 
dominant part of (|5.3p - is presented in Fig. [8l In the 
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(n+r) 



(J) 



(5.6) 



which leads to a universal exponent (|5.2p 

MJ) = 1. (5.7) 

In the deconfined regime < J < Jc we find a nonva- 
nishing singlet-triplet gap: 

^(n)_^(n) ^ ^(00) _ ^(00) ^ ^(„) ^g g^ 

with a finite-size correction 



(5.9) 



which follows from the difference r'^"+^) — of 
the largest eigenvalues t^"-+^\ jn the triplet and 

singlet sector and the fundamental coupling a*^") (|2.6p . 
which can be extracted from Fig. [4l The large n limit of 
the difference 



^(n+r) _ ^("+1) 



(5.10) 



turns out to be 2 for all J < J < Jc whereas the 
coupling 



a("+i)(J) = c^(J)-2-" 



(5.11) 



decreases with the system size iV = 4" as — 1= for all J 
< J < Jc- a{J) is shown in Fig. [9] 
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FIG. 9: n-dependence of a{J) = 2'" ■ a^"+'^> (J) - shown after 
n — 20 iterations. 



FIG. 8: The first derivative with respect to p. 

confined regime Jc < J < 1 all the curves approach a 



VI. THE STAGGERED MAGNETIZATION. 

Finally we want to present our results from the recur- 
sion formula [(8.3) in ref. dH)] 
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^(n+l) 



I |^(n+l) ^ 

= E 4"^'^<^f+'^r,,(7(")) (6.1) 

,i— 1 



for the staggered magnetization on an (n + l)-cluster. 
Note, that the renormahzation procedure only enters via 
the components a\"'^^\ i = 1, . . . , 7 on an (n + l)-cluster 
and the couphng 7^"^ 

The 7x7 matrix Tii ^{{j^"'^) is presented in Appendix 
C of (15). In Fig. [TD] we present the ratio as 
function of n and J for the case {do = 7, di = 9, d2 = 14) 
for large n; 



J=LO 
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J=0.4S3 
J=0.4S24 , 
J=0.4S23; 
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--- 0.25 ; 



JO ]2 14 16 18 20 



FIG. 10: The ratio as function of n and J for dimen- 

sions: do = 7, di — 9, d2 = 14. 

all the curves approach a common limit 

^ 0.7401.. for Jc < J < 1 (6.2) 
-> 0.25 for < J < 

VII. DISCUSSION AND PERSPECTIVES. 

We have studied in the 2D model with interacting 
plaquettes various observables like the scaled singlet- 
triplet gap (Fig. [3]) as function of n (i.e. sys- 
tem size 2"+^ X 2"+^) and the coupling parameter J in 
(jl.ip . We find spectacular differences in the confinement 
(Jc < J < 1) and deconfincmcnt (J < Jc) regime, which 
allows - for a given truncation scheme - for an extremely 
precise determination of the critical coupling Jc in all 
these quantities. This means, that the interaction ma- 
trices a1!'+^\ 5* = 0,1,2 (l^ - fTTU)) and thereby the 
renormahzation group equations [(6.2)-(6.6) in de- 
pend on J in an extremely sensitive way. The reason is a 
feedback between the scaled energy differences p*-"-*, k*^") 



- which enter in the diagonals of Ag 
eigenvalues cr("+i), t("+i', 



(«+i) 



- and the largest 



This feedback also leads to a dramatic change in the 
eigenstates |cr("+i)), |r("+i)), E.g. the square 

of the first component [ (|4.ip for i — I] (^cf^^^^^^ in 

the singlet eigenvector |cr'^"+^)) changes completely if we 
go from the confined (Jc < J < 1) to the deconfined 

(J < Jc) regime. In the deconfined regime (^<yi^^^^^ is 

almost one, which means, that the ground state factor- 

izes into 4 noninteracting singlets. In the confined phase 
2 

is very small. Therefore, the remaining com- 



(n+l) 



(n+1) 



i = 2, .., 7 contribute significantly to the 



ponents cr^ 
eigenstate {a'- 

These contributions are characterized by excitations 
of the cluster spins on the four plaquette system. Exci- 
tations of cluster spins are necessary to induce cluster- 
cluster interactions. The vanishing of the singlet-triplet 
gap - as it is observed in the confinement regime Jc < J < 
1 - is a consequence of the cluster-cluster interactions in- 
duced by cluster excitations (triplet and quintuplet). We 
have checked the dependence on the truncation of the in- 
teraction matrix by suppressing in Tables II, III, IV (of 
ref. (fish ) all states with one quintuplet plaquette. The 
dimensions of the interaction matrices reduce to 



{do,di,d2) = (5,3,4) 



(7.1) 



which of course worsens the renormahzation group ap- 
proach. This is signalled by a somewhat larger singlet- 
triplet gap. As a consequence the deconfined regime 
( J < Jc) is enlarged. 

If we look at the deconfinement parameter ^crj"^^^^ , 

Fig. [Tl]for the case (|7.ip . we observe a shift to a larger 
value of Jc'. 
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FIG. 11: The probability (^-["^^''^ evaluated for dimensions 
dLll). 



J^(5,3,4) = 0.5615.. ; Jc(7, 9, 14) = 0.4822.. . (7.2) 
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This value is close to the result of ref. ([loh obtained 
without quintuplet excitations. Therefore, the difference 
in the two values (|7.2p reflects the effect of rotational 
symmetric excited states on the 4 plaquette cluster with 
one quintuplet. We expect that further excited states 



with nq = 2,3,4 quintuplets will lead to changes in the 
values .has well. The Monte Carlo simulations of Janke 
et al. suggest, that the RG results should converge 
non-monotonously towards 0.549. 
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